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An equivalent-multishell approach for the approximate calculation of the characteristics of elec- 
tromagnetic waves propagating in almost circular (azimuthally symmetric), closely packed bundles 
of parallel, identical, and metallic carbon nanotubes (CNTs) yields results in reasonably good agree- 
ment with a many-body technique, for infinitely long bundles when the number of CNTs is moder- 
ately high. The slow-wave coefficients for azimunthally symmetric guided waves increase with the 
number of metallic CNTs in the bundle, tending for thick bundles to unity, which is characteristic of 
macroscopic metallic wires. The existence of an azimuthally nonsymmetric guided wave at low fre- 
quencies in a bundle of a large number of finite-length CNTs stands in contrast to the characteristics 
of guided-wave propagation in a single CNT. The equivalent-multishell approach yields the polar- 
izability scalar and the antenna efficiency of a bundle of finite-length CNTs in the long-wavelength 
regime over a wide frequency range spanning the terahertz and the near-infrared regimes. Edge 
effects give rise to geometric resonances in such bundles. The antenna efficiency of a CNT bundle 
at the first resonance can exceed that of a single CNT by four orders of magnitude, which is promis- 
ing for the design and development of CNT-bundle antennas and composite materials containing 
CNT-bundles as inclusions. 

PACS numbers: 42.70.-a, 73.25.+i, 77.84.Lf, 78.67.Ch 

I. INTRODUCTION 

Carbon nanotubes (CNTs) have been the subjects of intensive research in the area of nanotechnology for about 15 
years [l|, yet understanding of their fundamental physics is far from complete. One topic under intensive investigation 
is their electromagnetic response. Early theoretical studies on CNTs, modeled as infinitely long cylinders of a gas of 
electrons, showed the existence of gapless, low-frequency, plasmon branches The low-frequency plasmon excitation 
along the nanotube axis lead to the formation of a electromagnetic surface wave Q. This surface wave strongly 
influences the scattering [3] and radiation Q properties of a CNT in the terahertz regime. 

Guided-wave propagation should therefore occur in a CNT bundle containing 2 ^ 1000 parallel CNTs closely 
packed together. A guided wave could be formed by the plasmonic excitation in every CNT in the bundle and its 
characteristics would strongly depend on interactions between the CNTs. Two types of such interactions have been 
discussed in the literature: 

(i) The first type of interaction arises from the direct coupling of the electronic states in adjacent CNTs due to 
overlap of their electron wavefunctions j7|. The overlap is always very small for two reasons: (a) the contact area 
between two adjacent CNTs is small due to their geometric curvature, even when the two touch each other; and (b) 
the orbitals of the carbon atoms strongly overlap only in the plane of a graphene sheet, which leads to the known 
van der Waals form of the integraphene sheet interaction. The momentum mismatch between the Fermi points of 
neighboring CNTs suppresses the inter-CNT tunneling and leads to strong localization of electronic eigenstates on 
individual CNTs. Therefore, there is no significant change of low-energy band structure in the vicinity of the Fermi 
energy of a compositionally disordered metallic CNT bundle 1] as compared with band structure of a single CNT. 

(ii) The second type of interaction is an electrodynamic coupling in which Coulomb interactions in a CNT are 
modified by the dielectric screening induced by the adjacent CNTs. Such dielectric screening has a significant effect 
in ID structures (such as single-wall CNTs) since many-body interactions are intrinsically strong in ID geometry 
because much of the electric field of a charge on a CNT extends outside of the CNT. 

A bundle of parallel CNTs is generally considered as a 2D array of infinitely long CNTs [1, [^, [l3| , thereby allowing 
the determination of the dispersion properties of a low-frequency plasmon in a 2D periodic medium. A realistic bundle 
of almost circular cross-section has a cross-sectional diameter much less than both its length and the electromagnetic 
wavelength in free space. In order to determine the scattering and radiation properties of a realistic CNT bundle, one 
needs to consider both the finite diameter and the length of the bundle. 

Our aim in this paper is to analyze guided-wave propagation in an almost circular bundle of metallic CNTs all of 
which are either infinitely long or have finite length. We neglect modifications of the low-energy band structure of 
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CNTs, but not the electromagnetic coupling of CNTs in a bundle. Following Selpyan et al. [3|], we assume that, in 
the low-frequency regime below optical interband transitions, the conductivity of metalic CNTs is described by the 
electron-gas model, thereby leading to the CNT conductivity to follow the Drude model. 

The rest of this paper is organized as follows. In Sees. HIl and Hill a many-body technique and an equivalent-multishell 
approach, respectively, are applied to a bundle of infinitely long metallic CNTs to derive dispersion equations for 
guided-wave propagation on the bundle. In Sec. llVi scattering theory is applied to a bundle of finite-length CNTs. 
Sec. |V] contains numerical results for guided-wave parameters obtained from different approaches. 



II. GUIDED WAVES IN A BUNDLE OF INFINITELY LONG CNTS 



A. Many-body technique 

Let us examine the propagation of a guided wave in an isolated bundle of N metallic CNTs that are closely packed 
together and are of infinite length. The surrounding medium is free space (i.e., vacuum). The effective cross-sectional 
diameter of the bundle is much smaller than the wavelength in free space, and the transverse current in all CNTs in 
the bundle is neglected. An exp(— iwt) time-dependence is implicit, with i — \/— T, t denotes the time, and uj is the 
angular frequency. 

In the low-frequency regime, only intraband transitions of 7r-electrons with unchanged transverse quasi-momentum 
are allowed [lH . These transitions contribute to the axial conductivityjbut not to the transverse conductivity [13] , and 
excite azimuthally symmetric electric current densities in the CNTs [3, [l^ . At frequencies far away from interband 
resonances, in practical terms, azimuthally nonsymmetric electric current densities are not excited in CNTs because 
the relevant conductivities vanish. In detail, this peculiarity has been discussed elsewhere 0, with application 
to surface-plasmon modes. In the remainder of this paper, therefore we restrict ourselves to azimuthally symmetric 
electric current densities in the CNTs forming the bundle. 

The superposition of the fields, induced by the electric current densities on the surfaces of all CNTs together form 
a guided wave in the CNT bundle. The electric Hertz vector nm(r) = nm(r)ez created by the axial electric current 
density on the surface of the m*'' CNT (m G [1, -^]) in the bundle is governed by the Helmholtz equation 

(V2 + fc2)n„(r) =0, (1) 

where is the unit vector along the CNT axis (and therefore the axis of the bundle), k = lo/c, and c is speed of 
light in vacuum. Let the origin of the cylindrical coordinate system (p^™) , 0'^'") , z) be located at the point z = on 
the axis of the m*'' CNT. Since lim is a function only of p*^™' and z for an azimuthally symmetric electric current 
density, the nonzero components of the electromagnetic field in this cylindrical coordinate system are as follows: 



Uz2 



— ik 



(2) 

(3) 
(4) 



The electric Hertz potential H,„ must satisfy the effective boundary conditions 



dpi'") 



p('")=_R„+0 
k^Urr, 



H. 



dpi"" 



p('")=i?„„-0 
N 



(5) 



(6) 



where dm is the axial conductivity of the m*'* CNT in isolation. The scalar fields i?™"(z) and E'^^(z) are the z-directed 
components of the electric field on the surface of the m*'* CNT, produced by the externally impressed sources and the 
axial current density on the s*'' CNT, respectively. 
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Below the frequency regime of optical transitions, an expression for the axial conductivity (t„ 
quantum transport theory as 



is available via 



TT^hR^ (uj + i/t) ^ J dp^ dp^ 

IstBZ 



(7) 



where e is the electron charge, h is the normalized Planck constant, and Pz is the axial projection of the electron 
quasi-momentum. The integer n e [l,?7i] labels the 7r-electron energy bands, with m appearing in the dual index 
{rh,h) used to classify CNTs [llj. The time constant r of the electron mean-free-path is assumed to be equal to 
the inverse relaxation frequency. The abbreviation IstBZ restricts the variable Pz to the first Brillouin zone. The 
equilibrium Fermi distribution function 

FciPz, n) = -TT (8) 



1 -I- exp 



involves the temperature T and the Boltzmann constant kg. The electron energy SdPzi i^) for zigzag (m, 0) CNTs is 

El 



70 



/ nn \ 

1-1-4 cos (apz) cos ( + 4 



2 f'^'^ 
cos^ I 

V TO 



(9) 



where 70 ~ 2.7 eV is the overlapping integral a = 3b/2h, and b — 0.142 nm is the interatomic distance in 
graphene. An expression for SdPzTn) for armchair (to, to) CNTs is also available [ll[. 



B. Dispersion equation 

A solution of Eq. ^ that satisfies the boundary condition © is as follows: 

n,„(p('"),z) = 



\ /o(/«i?™)Xo(Kp('")), > R,n ■ 



Here is an amplitude, h is the guide wavenumber to be determined, and k — \/h'^ — k?^ while /o(-) and Ki^{-) are 
the modified Bessel functions of order 0. Equation ^ still has to be satisfied. 

The field E™'^{z) in Eq. (O may be found by applying Eqs. ^ and UHl) in the coordinate system (p*-*-*, z) to 
the s*'* CNT, s ^ TO, and then using the addition theorem for Ka{Kp^'^'>), p'^'*^ > Rg, to translate to the coordinate 
system (p^™), z) [B]. With p'-'^ lying on the surface of the m*^ CNT, as shown in Fig. [Ha), we get 

Erip^'lz) = -AsK^e''''Io{nRs)Ko{Kp('^) = 

00 

-AgK^e^'^^'IoinR,) ^ K,{nd,m)h{nRm)e''^'^\ (11) 

l=~oo 

where ds„i is the distance between the axes of the two CNTs and the angle (j)^"^'^ has been identified in Fig. [Ha). 

As the conductivities for azimuthally nonsymmetric modes {(. ^ 0) are assumed to be null-valued, only the I — Q 
term in Eqs. (|lip survives after the substitution that equation into Eq. ((S]). Thereafter, the substitution of Eqs. (flOl) 
and pT|) into Eq. ([5]) with E™^(z) = and the subsequent use of the Wronskian of modified Bessel functions [l^ 
lead to the following set of linear homogeneous equations with unknown As, s £ [1^ N]^ written in the matrix form as 

MA = 0. (12) 

Here, A is a column vector containing the N unknowns Ag, and the element Msm of the N x N matrix M is given by 

1^ ^ / Ko{Kdsm)lQ{KRm), S 7^ TO, 

\ Ko{KRs)~iuj/[4:TTRsCrsK^Io{KRs)], s = m. 

The set (|T2|) of linear equations has nontrivial solutions provided 

detM = 0. (13) 

The dispersion equation (|13p has N roots corresponding to N guided waves in the CNT bundle. The solution of 
Eq. [T51 allows us to obtain the slow- wave coefficients (3 = k/h oi guided waves in the CNT bundle. 
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FIG. 1: (a) Schematic of two CNTs in a bundle to identify various quantities appearing in Eq. (|f f [) . (b) A bundle of 
N = 55 closely packed CNTs that can be approximately taken to have azimuthal symmetry, and thus has an almost circular 
cross- section. 



III. EQUIVALENT-MULTISHELL APPROACH 



The many-body technique is cumbersome for a CNT bundle with large N. Furthermore, we are interested in guided 
waves with azimuthal symmetry, corresponding to the low-h roots of Eq. (|13p . as such waves are easily excited in 
bundles by uniform external fields, and are also relevant for bundles of finite-length CNTs in the long-wavelength 
regime (Sec. |TV|. 

An approximate but simpler method for \ow-h guided waves can be devised for an almost circular bundle of 
closely packed and identical CNTs as follows. Let Rt stand for an effective cross-sectional radius of the bundle. 
Since the electromagnetic field of the guided wave is azimuthally symmetric both inside and outside the bundle, 
we can model the bundle as a multishell structure comprising N concentric shells. Each shell in this equivalent- 
multishell structure is infinitesimally thin. Thus, the p*'' shell, p £ [1,^], has a cross-sectional radius Rp such that 
Rb = Rj^f > -R^v-r ••• ^1- "^^^ effective surface conductivity of the p*'' shell is given by dp = Ep/(27rJ?p), where 
is equal to the sum of the axial conductivities of all CNTs placed between the p*'* and the (p— 1)*'' shells. There 
is some latitude inherent in the procedure to select TV, Rp, and dp, p G [1,TV]. 

Let Tip{p,(t),z) be the Hertz potential everywhere entirely due to the p*'' shell, with the cylindrical coordinate 
system [p, (f), z) located at the point z = on the axis of the CNT bundle. The boundary conditions across the p*^ 
shell due to the axial current density 3p{z) — Jp(z)ez on the surface p = Rp are as follows: 



an. 



dp 



p=Rp+0 



dp 



ike ^ 



(14) 



p=flp+0 



p=B.p-0 



Here 



N r 



E 

9=1 



dz^ 



ng{Rp,^,z) + E°{2 



(15) 



(16) 



where i?^ is the z-directed component of the externally impressed electric field. 

With Ap^\ £ e [0, oo), representing its amplitude, an expression for Hp is as follows: 



Up{p,cl,,z)^4'h^'^^ 



Ki{nRp)h{Kp), p < Rp, 
Ii{KRp)Kt{Kp), p > Rp. 



(17) 



Substitution of with E'^ — and Eqs. (fTT]) into Eq. (|14p leads to a set of linear homogeneous equations with 
unknown p £ [I, N]. This set has nontrivial solutions that can be ascertained by solving dispersion equation 



det M = 



(18) 



for the determination of h. The element M„p of the N x N matrix M is to be computed as 



Kt{KRq)h{K.Rp), q<p, 
M„„ = { KiXnRp)h{K,Rq), q> p, 

Ke{nRq)Ie{KRq) - iuj /[iTTRqaqK^], q^p. 
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For almost circular bundles with iV = 55 CNTs, we compared the first three solutions of Eq. (|13p with those of Eq. 
(fT5)) for £ — 0, and obtained good agreement, as discussed in Sec. |Vl 

In contrast to an isolated CNT and even a bundle with relatively small number of closely packed CNTs, in the low- 
frequency regime an azimuthally nonsymmetric guided wave with £ can exist in a CNT bundle with a large number 
of CNTs. This guided wave is formed by the ensemble of azimuthally symmetric electric current densities excited in 
every CNT of the bundle. The local field of the azimuthally nonsymmetric wave quickly changes in the central part 
of the bundle. Since symmetric electric current densities are excited by the spatially homogeneous component of the 
local field, then CNTs in the central core of the bundle practically are not excited and their conductivity may be 
supposed to be zero. Therefore, in order to approximately describe an azimuthally asymmetric guided wave by Eqs. 
P?)) and the conductivity of the inner shells with radius Rp < lOiRo/ir should be assumed to be equal to zero, 
where Rq is the radius of every CNT in the bundle. It is expected that an azimuthally nonsymmetric guided wave 
can be easily excited in a bundle of parallel closely placed CNTs and contribute greatly to interbundle interactions. 



IV. GUIDED WAVES IN A BUNDLE OF FINITE-LENGTH CNTS 



In order to investigate the finite-length effects in CNT bundles, let us apply integral-equation methods developed 
for a single CNT [1, |^ and for a planar array of CNTs [H, [13] ■ Let a closely packed bundle of parallel and identical 
CNTs of length L be aligned parallel to the z axis of a Cartesian coordinate system. The bundle has an almost circular 
cross-section so the bundle radius Ri, can be prescribed; i?^ is assumed here to be much less than the wavelength A 
and length L. On exposure to an externally impressed field that is almost homogeneous over the bundle cross-section, 
an axial and azimuthally symmetric surface current density is excited in every CNT. Following Sec. IIIIl we replace 
the CNT bundle by N multishells, and prescribe the radius Rp as well as the effective axial conductivity dp of the 
p''' shell, p G [1, The surface current density, induced by the externally impressed field on the surface of the p"* 
shell, is denoted by Jp{z) = Jp{z)ez- K has to satisfy the edge conditions 



Jp(±i/2) = 0, 



(19) 



which express the absence of concentrated charges on the two edges of the bundle. 

The electric Hertz potential Ii{p,z) satisfies the Helmholz equation ([1]), the radiation condition [18], as well as the 
the boundary conditions (|14p and (fTS]) : the cylindrical coordinate system {p,(j>,z) located at the point z = on the 
axis of the CNT bundle. The potential n(p, z) is expressed in the form of a single-layer potential as 



n(p, z) = — / Jp{z')G{z - z' ,p,Rp)dz' . 



(20) 



-L/2 

where 

/^'^exp \ ik^ + R"^ — 2Rp cos tp + z'^ > 



is the free-space scalar Green function and 



+ _ 2Rp cos ip + z'^ 



(21) 



Jp{z) 



d'^n{R 



pi ' 



dz'^ 



k^U{Rp,z) + E°{z) 



(22) 



Setting p — Rp in Eq. and making use of Eq. (P^ . we obtain a system of N integral equations with respect 
to the unknown current densities as follows: 



*pW = E / y-^G{z-z\R„Rp) 



2ika, 



'-T- exp(ifc|z — z'\) 



xJq{z')dz', pe[l,N]. 



(23) 
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Here, 



1 



L/2 



2ifc(T„ 



E'^^{z')exp{ik\z- z'\)dz' 



-L/2 

+Cp exp(iA;z) + Dp exp{—ikz) 



(24) 



and Sqp is the Kronecker delta, whereas Cp and Dp are unknown constants to be determined from the edge conditions 
(fini) . ParentheticaUy, the system with a different Green function was appHed to a planar array of finite CNTs 

The integral on the right side of (j23[) can be numerically handled by a quadrature formula, thereby transforming the 
system (|23p into a matrix equation. The solution of the corresponding characteristic equation yields eigenfrequencies 
and eigcnmodcs of a finite coaxial cylinder as a high-Q microcavity. 

In the long- wavelength regime (A ^ L), the electromagnetic properties of the CNT bundle can be characterized by 
the polarizability scalar 



27ri 



L/2 



Jp{z)dz 



(25) 



As shown elsewhere [1, [E, @ , an isolated CNT can function as an antenna in the terahertz regime wherein the CNT 
has geometrical resonances of guided wave with slow-wave coefficient /3o at frequencies related to the CNT length L 
by the condition 



Lk^ nSReiPo), S=l,2,. 



(26) 



where /3o is the slow-wave coefficient for an isolated CNT. The antenna effect of an array of multiwall CNTs was 
experimentally found at a frequency satisfying the condition (pS)) with s — 1 and /3o ~ 1 [19, 20]. 

In this paper, we are interested in the antenna efficiency of a CNT bundle at the first antenna resonance: s = 1 in 
Eq. The antenna efficiency is defined as the ratio 



Pr 



Pt+Pr 



where 



Pr = 



sm 



L/2 



N 

*'=^™"^^J?pJp(z)dz 
-L/2 „=i 



de 



(27) 



(28) 



is the radiated power and 



Pf 



/ W n rL/2 

TTRe y ^ / \Jp{z)\^dz 

\t^l J -L/2 



(29) 



is the power lost to ohmic dissipation. 

In the long-wavelength regime, the antenna efficiency 770 of an isolated CNT is equivalent to that of a dipole antenna 
of length L and resistance per unit length Rdip 0; i-S., 



L/X 



L/A-t-3ci?dipA/(87r2 



(30) 



where Rdip = Re[l/(27ri?oO'o)] for an isolated CNT of cross-sectional radius Rq and surface conductivity (Tq- The high 
value of Rdip p3| and the small value of L/ A at the first antenna resonance lead to a very small antenna efficiency of 
an isolated single CNT, i.e., 770 « 10"^ - 10"^ [li]. 

However, the situation is more optimistic for a bundle of N metallic CNTs. Such a CNT bundle can be considered 
as a composite antenna containing N in-phase parallel dipole antennas, so that Pr ^ N"^ but Pt ^ N from Eqs. ((28|) 
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and (|29p. Then the antenna efficiency of the CNT bundle at the frequency of first geometrical resonance of the axially 
symmetric guided wave with the highest slow-wave coefficient (3 is 

^^oR^(/^) (31) 

' 77o[AfRe(/3)-Re(/3o)] + Re(/3o)' ^ ^ 

where Tyo is the antenna efficiency of an isolated CNT at the first antenna resonance at the same frequency. 

Because of electromagnetic coupling of the CNTs in a bundle with 1 CNTs, the inequality Re(/9) ^ Re(/3o) 

holds true, as shown in Sec. |Vl therefore, r/ 3> 770- 



V. CHARACTERISTICS OF GUIDED WAVES 



Calculations were performed for almost circular bundles made of N parallel, identical, single-wall, metallic, zigzag 
(21,0) CNTs. The relaxation time was taken as r = 10~^^ s. For convenience, the CNTs were assumed to be 
close-packed on a 2D triangular lattice with intertube spacing 3.4 A as shown in Fig. [ijb) for N — 55. 



A. Guided waves in bundles of infinitely long CNTs 




FIG. 2: Radial dependencies of the magnitude of the z-directed component of the electric field (in arbitrary units) of the guided 
waves in a CNT bundle of = 55 parallel, identical, infinitely long, single-wall, metallic, zigzag (21,0) CNTs arranged on a 
triangular lattice. The points correspond to the solutions of Eq. (|13p . and the lines to the solutions of Eq. (|18|) . 

Let us begin with numerical results for bundles of infinitely long CNTs. We considered only three roots each of the 
dispersion equations ([13]) and (fT8| . These roots — labeled /ii, /12 and /13 — were the ones with the smallest real parts 
(Re(/i3) > Re(/i2) > Re(ft-i)), and correspond to azimuthally symmetric guided waves identified as SW\, SW2 and 
SWi, respectively. Only these guided waves mostly infiuence the scattering properties of finite-length CNT bundles, 
as discussed in Sec. IV Bl We also considered the properties of an azimuthally nonsymmetric guided wave, identified 
as NSWl, which emerges from the solution of Eq. (fT8)) for i — \ and the real part of whose wavenumber h is the 
smallest possible. 

The radial dependencies of the magnitude of the z-directed component of the electric field inside the chosen CNT 
bundle for guided waves SWl^ SW2^ SWi and NSWl are shown in Fig. O The lines in this figure were obtained 
from the effective multishell approach of Sec. Illli whereas the points were obtained from the many-body technique of 
Sec. Ill Bl We conclude that the two approaches yielded reasonably close results for the chosen bundle. Such a good 
coincidence of the results from both approaches was observed only for bundles with closely packed CNTs. For rarified 
bundles (where the smallest inter-CNT distance exceeds the CNT radius by a factor of 10 or more), the results of 
both approaches coincide only for the guided wave SWl. 

As one can be deduced from Fig. ^ the electric field inside the bundle is distributed over the entire cross-section. 
Outside the bundle, the radial distribution of the electric field is governed by the decreasing modified Bessel function 
Ke{\/h'^ — k'^p). Thus the electric field is highly localized to the CNT bundle. 
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FIG. 3; Radial dependence of the axial component of the time-averaged Pointing vector (in arbitrary units) of the guided 
waves in the same CNT bundle as in Fig. (2] 




FIG. 4: Dependencies of Re(/3) of the guided wave SWl on (a) the number A'^ of CNTs in the bundle and (b) the bundle radius 
Rb- 

Furthermore, calculations show that the radial field distribution for SW2 and SW3 can be described adequately by 
Bessel functions of the first kind Jo(k2P) and Jo(k3p), respectively, where K2 and K3 > K2 are the non-zero minimal 
roots of the equation Ji{nRb) = 0. 

All other results — presented in Fig. [3][7] — were obtained using the effective multishell approach. Figure [3] shows 
the radial distributions of the axial component of the time-averaged Pointing vector Sz = c\H ^\'^ / for three 
azimuthally symmetric guided waves {SWl, SW2, and SWi) inside the CNT bundle {p/Rb < 1) and in the vicinity 
of its surface {p/Rb > !)• As shown in Fig. [31 the axial component of the time-averaged Pointing vector of 5*1^1 
is maximum near the surface of the bundle, and a large part of the energy of SWl leaks outside the bundle. In 
contrast, the power densities of SW2 and SW3 are mostly concentrated inside the bundle. Thus, the electromagnetic 
energy and volume electric current density of the guided waves in a CNT bundle are distributed similarly to those of 
eigenwaves propagating in a macroscopic, lossy, infinitely long wire [23|. Accordingly, guided waves propagating in a 
CNT bundle with a large number of CNTs cannot be called surface waves. 

Figure m^a) contains plots of the real part of the slow- wave coefficient /3 of the guided wave SWl in the terahertz 
{v = uj/2tt = 2.5 THz) and infrared {i/ = 27 THz) regimes for different numbers of CNTs in the bundle. The coefficient 
Re(/3) increases as much as 26 times with N increasing up to 900. A large N means that the bundle is thick (with 
Rb > 25 nm) and its slow-wave coefficient tends to that of a macroscopic metallic wire. The dependence of Re(/3) on 
the bundle radius is linear up to Rt = 25 nm, as may be deduced from Fig. [Hb). 

A comparison of Figs. [H^a) and Fig. [5] reveals that Re(/3) of SWl is about three, five and nine times more than 
that of NSWl, SW2 and SW3, respectively. We also found that the solution of the dispersion equation (flS)) for 
NSWl does not change in a frequency regime wherein the arguments of Bessel functions (•) and Ki{-) in Eq. (fTT]) 
are much smaller than unity, i.e. kRi, << 1. This condition holds for realistic CNT bundles over a wide frequency 
range from the terahertz to the near-infrared regimes f^. The dependencies of Re(/3) on Rb are linear for SW2, SW3 
and NSWl are also linear up to Rb = 25 nm (not shown in this paper). 
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FIG. 5: Dependence of Re(/3) on N, for the guided waves SW2, SW2 and NSWl at v = uj /2tv = 27 THz. 



Although not substantiated here by a graph, the value of — Im(/3)/Re(/3) for the considered guided waves does not 
depend on N. This implies that the electromagnetic coupling of the CNTs does not influence the attenuation of the 
guided waves in the bundle. The ratio — Im(/3)/Re(/3) is approximately equal to 0.32 and 0.034 aX i' — 2.5 THz and 
1^ = 27 THz, respectively, for SWl. 




FIG. 6: Frequency dependence of (a) lm{a) and (b) Re(Qf) of a bundle of CNTs of lengtli L — 500 nm. Two values of A'' were 
used: 55 and 800. The labels 1, 2 and 3 denote the geometric resonance of the guided waves SWl, SW2 and SW3, respectively. 
The subscripts on the labels denote the number s of the geometric resonance in Eq. (|26p . 



B. Guided waves in bundles of finite-length CNTs 

Let us now move on to almost circular bundles of finite-length CNTs. Fig. [S]demonstrates the frequency dependence 
of the polarizability scalar a of a bundle of CNTs for two different values N. The labels 1, 2 and 3 in this figure 
denote the geometric resonance of the guided waves SWl, SW2 and SW3, respectively. Clearly, the polarizability 
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resonances in this figure occur at frequencies satisfying the condition (|26|) . The location of the first resonance (s = 1) 
of all three guided waves on the frequency axis depends on N. The first geometrical resonance of the SWl is the 
strongest of the three, and it shifted from the terahertz regime to the mid-infrared regime as N was changed from 1 to 
800. The geometrical resonances of SW2 and SW3 occur at lower frequencies; moreover, they even vanish for small 
A^, because of the strong attenuation of the guided wave at low frequencies (where the condition Im(/i)/Re(/i) > 1 
holds). 




FIG. 7: Dependence of the antenna efficiency -q on the number A'^ of CNTs ai u — 2.5 THz (terahertz regime) and f — 27 THz 
(infrared regime). The length L was modified to ensure that the condition (|26|l is satisfied with s = 1, for fixed u. 

Similar conclusions could be made for wavelength-dependent resonant effects in the optical regime from measure- 
ments on a parallel array of multiwall CNTs [l^. Thus, by varying the bundle radius and the CNT length, one can 
tune the resonance properties of a CNT bundle functioning as an antenna from the terahertz to the visible regimes. 
A similar conclusion could be drawn for composite materials containing CNT inclusions [1, [2^ . 

The dependence of the antenna efficiency rj on the number N of CNTs in a bundle at two different frequencies 
1/ — ioj^-K is illustrated in Fig. [T] The length L was modified to ensure that a chosen value of v always corresponds 
to the first geometrical resonance of SWl (i.e., the condition ([^5]) is satisfied with s = 1). Evidently from Fig. [71 t] 
increases with N and tends to unity for thick bundles {N > 800); indeed, rj — 0.24 when N — 900 and v = TI THz. 



VI. CONCLUDING REMARKS 



To conclude, an equivalent-multishell approach was proposed for the approximate calculation of the characteristics 
of electromagnetic guided waves on almost circular, closely packed bundles of parallel, identical, and metallic carbon 
nanotubes. The CNTs can be either infinitely long or of finite length. The dispersion characteristics of the guided 
waves with the smallest retardation (i.e., lowest values of Re(/i)) where studied for bundles of infinitely long CNTs. 
The slow- wave coefficients for azimunthally symmetric guided waves were found to increase with the number of metallic 
CNTs in the bundle, tending for thick bundles to unity, which is characteristic of macroscopic metallic wires. The 
existence of an azimuthally nonsymmetric guided wave at low frequencies in a bundle of a large number of finite- length 
CNTs was demonstrated, in contrast to the characteristics of guided-wave propagation in a single CNT. 

The polarizability scalar and the antenna efficiency of a bundle of finite-length CNT in the long-wavelength regime 
were calculated over a wide frequency range spanning the terahertz and the near- infrared regimes. The resonances of 
different guided waves in CNT bundles caused by edge effects (geometrical resonances) were identified. The antenna 
efficiency of a CNT bundle at the first antenna resonance can greatly exceed that of a single CNT. Thus, the analysis 
carried out in this paper forms a basis for the design and development of CNT-bundle antennas and composite 
materials !24] containing CNT-bundlcs as inclusions. 
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